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ABSTRACT 

We have studied the effects of Compton cooling on cooling flow by performing numer- 
ical hydro dynamic calculations of the time evolution of hot gas in clusters of galaxies 
with luminous quasars. We assumed various temperatures for the hot intracluster gas. 
We have shown that the Compton cooling due to very luminous quasar is effective in 
inducing cooling flow, before radiative cooling flow is realized. The mass flux due to 
the cooling flow increases with time, as the Compton cooled region expands. However, 
the mass flux of the Compton cooling flow is not large, less than 1 Mq yr _1 in our 
model, since Compton cooled region is limited in an inner galactic region around a 
quasar. Even though the quasar active phase ceased, the cooling flow will continue 
for at least 10 9 yr. The accreted mass is enough to explain X-ray absorption lines 
in high red shift quasars, if the Compton cooled gas is compressed by high pressure 
intracluster gas. 
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1 INTRODUCTION 



Radio-loud quasars have some characteristic features differ- 
ent from radio-quiet quasars. Radio- loud quasars are ob- 
served to be strong UV sources (Luv ~ 10 46 — 10 47 erg s 
_1 ) (e.g., McDowell et al. 1989), and host galaxies of radio- 
loud quasars are found preferentially in high density envi- 
ronments (e.g., Smith and Heckman 1990) and are often 
central elliptical galaxies in clusters (Malkan 1984). Cool- 
ing flows sometimes appear arou nd th ese quasars (Fabian 
1994). For examples, Elvis et al. (1994) have shown strong 
X-ray absorption lines in high redshift quasars. They have 
suggested that the origin of absorbing gas is a cooling flow 
around quasar and have estimated that the column density 
of gas that absorbs the X-ray is up to 10 23 cm -3 . Usually 
cooling flows in clusters of galaxies are proposed to be ac- 
cretion flows of intracluster medium induced by radiative 
cooling (Sarazin 1986; Fabian et al. 1987; Forman 1988). 
Since radiative cooling time may be longer than the cos- 
mic age at such a high z, another cooling process is needed 
for the cooling flow around quasar. Authors suggested that 
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cooling flows are induced by Compton cooling due to strong 
UV photons of quasars. 

For the majority of quasars the bulk of the emitted elec- 
tromagnetic luminosity is observed in the optical/UV/soft- 
X-ray spectral domain, the so-called big blue bump (Shields 
1978; Malkan 1983). The origin of the big blue bump is accre- 
tion through a Keplerian disk or a non-Keplerian torus onto 
a supe rmassive black hole (Rees 1984) . Fabian and Crawford 
( 1992| ) (hereafter F.C.) have proposed that, since Compton 
scattering between strong UV photons and hot electrons can 
dominate cooling process for the hot gas in central region 
of a cluster, accretion by Compton cooling can supply the 
gas to a central massive black hole, and have shown that 
the accretion rate by Compton cooling flow can be larger 
than that of Bondi-accretion flow. We summarise this point 
according to F.C. as follows. 

The Bondi accretion rate, Mb, may be estimated as 



Mb 



-(GM) p, 



(1) 



where p, c s , M, a, are the gas density, sound speed, the cen- 
tral black hole mass and a factor depending on the adiabatic 
index of the gas, respectively. We estimate the typical value 
of M B in equation <JlJ> as, M b ~ 2.0 x 10~ 3 Mq yr" 1 , where 
we take p = 6 x 10~ 26 g cm -3 , temperature T = 10 7 K, 
M = 10 9 Mg Mq and a = 1. At a radius R, the time-scale 
by Compton cooling, r c , is 
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r c « 10 13 R\ L^ 1 s, (2) 

where the radius 7? = 100 R2 pc and the luminosity L = 10 47 
L47 erg s . While the time-scale by bremsstrahlung cooling, 

Trad, is 

r rad Ri6x 10 13 T 3/2 Ps 1 s, (3) 

where the temperature T = 10 7 T7 K and the pressure P = 
nT = 10 s Ps cm -3 K. Therefore within a radius 



Rc « 250 L\(? T 3/4 P 8 ~ 1/2 



pc, 



(4) 



Compton cooling dominates bremsstrahlung cooling in the 
highly ionized conditions, so that we can expect stable gas 
inflow since the cooling time does not depend on gas density. 
Consequently, all the gas at R c flows into the central black 
hole, the accretion rate in this case is 



M c 



An Rip 

'fad 

1.45 L% 2 T 7 1/4 ( 



6 x 10" 26 g cm~ 3 ' 



a/a 



M, 



yr 



"(5) 



that is much larger than the Bondi accretion rate. The de- 
pendence of M c on L, i.e. M c <x L 3//2 and the linear de- 
pendence of L on M, means that, if L > L c , the accretion 
induced by Compton cooling increases L further, where L c 
is the quasar luminosity at the accretion rate of equation 
(^). On the other hand, if L < L c then M c is reduced and 
the luminosity continues to decline to Lb that is denned by 
the luminosity of a quasar at the mass accretion rate of Mb . 

When Compton cooling time is shorter than radiative 
cooling time in the inner central region, Compton cooling 
can induce the gas inflow before the radiation cooling in- 
duces the gas inflow. However, when the life time of quasar 
is comparable to Compton cooling time of the gas within R c , 
gas flow can not be steady state during the duration time of 
the quasar. Therefore we should study numerically the time 
evolution of the gas inflow induced by Compton cooling. 
Moreover we expect another possibility that by the density 
concentration due to Compton cooling radiative cooling time 
would be shorter than Compton cooling time, so that the hot 
gas would be cool radiatively at the central region and that 
the density concentration could go on. We examine whether 
the mass accretion rate increases by the effect of Compton 
cooling and how Compton cooling affects the evolution of 
the intracluster gas by using hydrodynamic scheme. We fix 
the quasar luminosity in our model for simplicity through 
our all simulations. 

In §2, we describe our model of the luminous quasar 
with the cooling flow in the central galaxy in a cluster. Our 
model incorporates the elliptical galaxy that includes stars, 
dark halos and cooling flow. In §3, we describe the results 
of our model calculations. In §4, we discuss the results and 
summarize our conclusions. 



2 A MODEL FOR COOLING FLOW WITH 
LUMINOUS QUASAR IN A CENTRAL 
CLUSTER GALAXY 

The central concept of our model is the evolution of density, 
temperature and accretion rate of the cooling flow in lumi- 
nous elliptical host galaxy with a quasar. We assume that 



the age of quasar is 10 yr. By treating the quasar as a source 
of UV photons, we consider Compton cooling process. 

Since we concentrate on the evolution of the intraclus- 
ter gas near the central galaxy, we calculate gas dynamics 
near the central galaxy. We assume that the cooling flow is 
spherically symmetric for simplicity. 



2.1 Mass distribution of the stellar system and 
dark halo 

The stars in the central galaxy are assumed to be distributed 
as the modified King profile with the density distribution 
given by 



p*(r) = p c [l + (r/r c ) 2 ]- 3 / 2 , 



(6) 



where r is the radial distance and r c = 500 pc is the core 
radius of the stellar distribution and p c = 3.0 x 10 -21 g 
cm -3 . Observations indicate that elliptical galaxies do not 
have flat cores, but rather cusp profiles (Hernquist 1990). 
However, the mass inside cusp is not so large, if we compare 
the mass within the core region. As we will show in Section 3, 
the cooling flow flux is mainly determined by the difference 
of pressures between at the inside and at the outside of the 
cooled region. Therefore, the choice of mass model with or 
without cusp is not so serious for our model. 

We have assumed that the dark halo has a total mass 
which increases linearly with radial distance as suggested 
for early- type galaxies by X-ray mass determinations (e.g., 
Forman et al. 1985). The particular dark matter density 
distribution we have adopted is 



Pd(r) = Pdc [ 1 + (r/r d ) 2 ]' 



(7) 



where r<j = 10 kpc is the core radius of the dark halo distri- 
bution and pdc = 6.26 x 10~ 25 g cm -3 . 

The stellar and dark matter mass distributions, M,(r) 
and Mu(r), can be determined by integrating the respective 
expressions for the densities. Carrying out the integrations, 
the mass distributions are found to be 



M»(r) =4tt p c r 3 [ ln(x + \/ x 2 + 1) - x/^fx 2 + 1 ] 



and 



Md (r) = 47r p dc rl (y - tan 1 y) , 



(8) 



(9) 



where x = r/r c and y = r/rd- The corresponding virial 
temperature of the model galaxy at r c is 2 x 10 7 K. 

We consider a supermassive black hole (hereafter 
SMBH) at the centre of galaxy. We assume that the mass 
of SMBH is M BH = 10 9 M . The SMBH mass dominates 
within the radius 



Mbh a 



^ = 172 < ire;) 173 pc 



At such a radius the virial temperature Tbh due to the 
SMBH is 



2.2 The initial conditions 

Given the mass distribution of the stars and dark matter, 
the initial density distribution of the hot gas can be derived 
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by assuming that it is in hydrostatic equilibrium. This as- 
sumption is valid when flow velocity is small. This condition 
can be written 



dP(r) 
dr 



-p(r 



GM(r) 



(10) 



where P(r) is the gas pressure, p(r) is the gas density distri- 
bution, and G is the gravitational constant. The total mass 
M(r) is given by 



M(r) = M«(r) + M D (r) + M B 



(11) 



The gas density distribution for an isothermal temperature 
is given by 

,(itriH , 



p(r)=p s0 exp[-(0(r)-0(O))- 



kT 



(12) 



where k is the Boltzmann constant, To is the initial isother- 
mal temperature, and 4>{r) is the gravitational potential by 
the stellar and dark matter mass distributions. We take 
p g o = 6 x 10 -26 g cm -3 , which maybe a typical value 
of interstellar matter in the host galaxy. We assume var- 
ious value of To as shown in Table 1, since we want to 
study Compton cooling effects on intracluster gas in var- 
ious types of cluster of galaxies. We redefine R c by using 
equation (2) and r c = 10 8 yr. By the new definition, we get 
R c = 540 pc for L uv = 10 46 erg s" 1 and R c = 1700 pc for 
Luv = 10 47 erg s" 1 . 

Since we assume that flow velocity is small except for 
the vicinity of a SMBH, the initial velocity distribution in 
the central region is given by 



v(r) — 



n e (r) ni(r)A r a r 2 



(13) 



p{r)GM{r) 

where n c (r) is the electron density, rii(r) is the ion density, 
Ao Tq is the emissivity. We use the steady state energy 
equation. In the outer region, the velocity distribution is 
given by 



v(r) 



■ ' ^(rinner) 



(14) 



r 2 p(r) 

We assume that rj nncr is r<j. 

We assume that cooling function A(T) = Ao T a . We 
t ake A (T) to be that given by Raymond, Cox and Smith 
(1976) for an optically thin, fully ionized gas with the cosmic 
abundance for T ^ 10 4 K and zero for T < 10 4 K. A(T) is 
given cgs unit by 



A(T) = 



2.2 x 10 

5.3 x 10 
6.2 x 10 
1.0 x 10 



-27 



T 

17 rp 
19 rp 
24 y0.55 



0.5 
0.87 
0.6 



for T > 4 x 10 7 
for 4 x 10 7 > T > 2 x 10 6 
for 2 x 10 6 > T > 10 5 
for 10 5 > T > 10 4 . 



dp 19,2 \ 



0, 



dpv 18,22, dP 



GM 



(16) 



(17) 



where 

p _ pkT 
(linn 

is the gas pressure and 

U = -pv H - 

2 7—1 

is the sum of kinetic energy and thermal energy per unit 
volume, and 



(19) 



(20) 



A C (T) = A c o L T 



(21) 



is the coefficient related with the Compton cooling and L is 
the UV luminosity of the quasar. A c o is the constant, 

1.1 a T k 



A c0 



(22) 



(15) 



2.3 The basic equations 

We calculate the time evolution of spherically symmetric 
X-ray emitting hot gas models. We assume that the gas is 
ideal. The hydrodynamic equations for a hot gas are 



2.1 (i m e c 2 ' 

where ctt is the cross section for the Thomson scattering. 
We take the 7 = 5/3. Other physical variables have usual 
meanings. 

2.4 Timescale 

We have performed numerical calculations in the variety of 
cases. In Table ??, we show the numerical model parameters 
and Bondi accretion rate Mb to compare with our numerical 
results. We take a = 1 and r c = 10 8 yr. In Fig. |l|, we show 
the cooling time (r c and r ra a) in the variety of temperature 
and luminosity. It is interesting that Compton cooling time 
depends on the radius and luminosity and does not depend 
on the gas density and temperature. 



2.5 The simulation 

We assume six cases of temperature of intracluster medium. 
Instead of the gravitational confinement of the hot intra- 
cluster gas in the cluster of galaxies, we assume the outer 
boundary confines this hot gas. We use the flux-split scheme 
with second-order accuracy in space and the first-order ac- 
curacy in time (van Leer 1982, van Albada, van Leer and 
Roberts 1982 ) . This scheme is the one of the upwind scheme 
for the Euler equation. We calculate the region from r = 50 
pc to 50 kpc. The calculated region covers the black hole 
mass dominated region as shown in Section 2.1. We use a 
mesh with 4996 elements. 

We take the boundary conditions as follows. At the in- 
ner and outer boundaries, the spatial gradient of the velocity 
is zero. At the inner boundaries, we assume that gas mass 
flux and total energy flux are same those in the cell next to 
the inner boundary. By this inner boundary condition, we 
obtain smooth accretion flow. At the outer boundaries,the 
density and energy are fixed in time. 

We have calculated the evolution of the cooling flow for 
the both cases of the existence of Compton cooling and no 
Compton cooling (only radiative cooling) for t — 1 x 10 9 yr 
from the initial situation. 



8U 1 d , 2 ,rr m , GM . ^ a pA c (T) . . 

o~{ r (U+P )«} = ~pv— 9 — n c mA T a - H c ; 1 , 18 
at r 1 or r 1 r 2 



3 NUMERICAL RESULTS 

We have shown that Compton cooling is dominant cooling 
process within R c which is defined by equation (4) and that 
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there is the region where Compton cooling time is shorter 
than the age of a quasar as shown in Fig. 1. In this region, 
we can expect a cooling flow induced by Compton cooling 
within the lifetime of a quasar. From our numerical results, 
we found that Compton cooling induces a cooling flow, even 
if the temperature of intracluster gas is as high as 10 s K. 
However, the mass flux is smaller than that estimated by 
equation (5), since the lifetime of a quasar is too short to 
develop the cooling flow structure. We found the case in 
which the cooling flow induced by Compton cooling con- 
tinues after the lifetime of quasar in the lower temperature 
intracluster gas case. 

In Fig. we show the results of 7.d747 model as the 
strong Compton cooling flow case. These figures show the 
evolution of the density, the velocity, the temperature and 
the accretion rate of the gas in 7.d747 model . In this model, 
we assume the luminous quasar with Luv = 10 47 erg s _1 
and the temperature of intracluster gas, 7 x 10 7 K. In this 
case, since Compton cooling is very effective, the gas in the 
inner region of galaxy cools by Compton cooling within the 
lifetime of a quasar and the cooled gas inflow is realized as 
shown in Fig. 2(b). We should note that radiative cooling 
time of the hot gas initially filled in the host galaxy is longer 
than the calculation time and we do not find any cooling 
flow in 7.d7rad model in which we do not consider Compton 
cooling. The inflow velocity attains 40 km s _1 at r c at 10 8 
yr. The mass flux is as large as 0.45 Mq yr -1 in this stage. 
This is smaller than that estimated by equation (5). After 
the lifetime of a quasar, we stop Compton cooling, that is, 
we calculate the energy equation without Compton cooling. 
After then, the cooled gas flows into the central region, and 
the hot gas occupies the whole region. Finally, the cooling 
flow disappears. 

Next, we show the weak Compton-cooling-flow case in 
Fig. 3. In Fig. 3, we show 7.d746 model. Fig. 3 shows that the 
temperature of gas within 500 pc decreases due to Comp- 
ton cooling. As a result, the gas density increases by the 
compression due to the high pressure of the outer hot gas. 
However, the flow velocity and the mass flux are very small 
in the inner region, as shown in Fig. 3(b) and Fig. 3(d). The 
hot gas in the inner region is in nearly hydrostatic equilib- 
rium. The small mass flux would not affect the evolution of 
a quasar and the host galaxy. Since the gas cooled by Comp- 
ton cooling is not thermally unstable as shown by equation 
(2) and the gas temperature in this region remains rather 
high temperature in comparison with the virial temperature 
of the host galaxy, the cooling flow is not realized. This is 
the reason why gas is nearly hydrostatic state. 

Next, we show the Compton-cooling-induced-cooling 
flow case in Fig. 4. Fig. 4 shows the numerical results of 
5d747 model and that the gas in the inner region rapidly 
cools by Compton cooling, and attains less than 10 7 K 
within 10 s yr in the central region of the host galaxy. In 
this case, if we do not consider Compton cooling, radiative 
cooling flow does not appear as shown in Fig. 5. Since the 
temperature of gas near r c becomes much less than the virial 
temperature of the host galaxy, 2 x 10 7 K, the gas inflow is 
induced and the mass flux attains to 0.4 Mq yr -1 . After 
10 s yr -1 , we stop Compton cooling, and the inflow velocity 
decreases. However, radiative cooling flow continues to 10 9 
yr. The mass flux is 0.2 Mq yr -1 at 10 9 yr. We can con- 
clude that the cooling flow after the lifetime of a quasar in 



5d747 model is initially induced by Compton cooling and is 
maintained by radiative cooling. 

Next, We show the numerical results of 5d746 model 
in Fig. 6. In 5d746 model, we assume one-tenth of a quasar 
luminosity of 5d747 model, 10 4 6 erg s _1 . In this case, the 
gas in the inner region of the host galaxy cools by Compton 
cooling and the density in the cooled gas region increases. 
However, the accretion rate is as small as 0.01 Mq yr -1 , 
since the Compton cooling rate is smaller than 5d747 model. 
After Compton cooling stopped, the gas density and temper- 
ature return to the initial state. Although the gas cools by 
Compton cooling, the cooling flow disappears after Comp- 
ton cooling stops. From these results, we can interpret that 
the strong Compton cooling induces radiative cooling flow 
in 5d747 model. 



4 DISCUSSION 

We have shown that, if a quasar has very strong UV source 
(Luv = 10 47 erg s^ 1 ), the cooling flow can be induced by 
the quasar during the quasar active phase. In this case, intr- 
acluster gas is accumulated into the host galaxy. The total 
mass accreted into the host galaxy is not so large but enough 
mass to maintain the activity of an active galactic nuclei. On 
the other hand, for the typical lifetime of a quasar, 10 8 yr, 
a quasar activity ceases before Compton cooling becomes 
effective to induce cooling flow for the high temperature in- 
tracluster gas case (e.g., initial gas temperature of ~ 10 8 K), 
if a quasar does not have a strong UV source (Luv = 10 46 
erg s" 1 ). 

Due to Compton cooling, the gas temperature de- 
creases, and the gas density increases by the compression 
of high pressure of outer gas. The gas flows in this stage 
corresponds to the pressure driven accretion (Sarazin 1986). 
If the temperature of cooled gas becomes lower than the 
virial temperature of the host galaxy, the inflow induced by 
Compton cooling occurs. As shown in Table 2, the mass 
flux at 10 s yr can be approximated by using the equation 
(5), except for replacement of r ra( j by 10 s yr and R c by 
the distance where Compton cooling time equals to 10 8 yr. 
Then, equation (5) should be changed as M c ~ 0.5 L 3 J 7 2 . 
It is possible, that the feedback on central accretion due 
to the increased quasar luminosity, as Compton cooling be- 
come effective, would increase the mass flow from large radii 
to large enough values to account for central accretion rate, 
if the quasar UV luminosity reaches 10 48 erg s _1 . 

Compton cooling affects the halo gas in the host galaxy 
and increases its gas density. As a result, after a quasar ac- 
tivity ceases and Compton cooling become ineffective, there 
are the cases that radiative cooling flow continues, although 
the mass flux is smaller. 

Compton cooling is not thermally unstable, since 
Compton cooling does not depend on the gas density. Al- 
though Compton cooling time is much smaller than the sim- 
ulation time, we have gotten the results that the tempera- 
ture of hot gas does not attain very low temperature in the 
inner region in the high temperature case. The reason is as 
follows: In 7d746 model, the flow time at 10 8 yr is ~ 10 7 yr 
at 100 pc, since v ~ 10 km s . There is not enough time 
to cool further for Compton cooled gas in the calculated re- 
gion, since Compton cooling time in this region is as large as 
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10 7 yr. On the other hand, in 5d747 model, the flow time at 
10 s yr is as large as 4 x 10 7 yr at r ~ 1000 pc, since v ~ 25 
km s _1 . Since Compton cooling time is as large as 10 7 yr 
at r ~ 1000 pc, Compton cooling is effective to reduce the 
temperature of intracluster gas. 

Elvis et al. (1994) showed the evidence of X-ray absorp- 
tion in high redshift quasars. They found that these quasars 
are very luminous X-ray sources up to 10 48 erg s" 1 . They 
suggested that a cooling flow around the quasar accounts 
for the X-ray absorption. If these quasars have comparable 
magnitude of UV flux with X-ray, Compton cooling would 
be very effective in these quasars. As discussed above, since 
M c oc I/uv 3 ^ 2 , the mass accretion rate becomes large. This 
cooled gas in Compton cooling flow would correspond to X- 
ray absorption. Wc compare column density of our Compton 
cooling flow models with the column density observed in the 
high redshift quasars. As shown in Section 3, the gas den- 
sity distribution may be approximated by p — poir/ro) -1 , 
where po ~ 2 x 10~ 25 g cm -3 and ro ~ 1000 pc. The column 
density N is estimated as 



N 



10 21 ( 



Po 



2 x 10~ 25 g cm 



-)( T 



r 



i In 



(r 



kpC (fin/1 pc) 



/I kpc) _ 3 
cm . 



This column density is smaller than the observed value. The 
compression of cooled gas is expected, since the pressure of 
central gas is very high. It is possible to explain the high col- 
umn density at quasar observed in high redshift by Compton 
cooling flow by this compression. 

In our numerical results, the cooled gas flows though 
the inner boundary. It is not clear that how the gas ac- 
cretes onto the central black hole or an accretion disk around 
the SMBH. Then we hold our model luminosity fixed. The 
cooling of inflow from large radii to accretion at small radii 
is beyond the scope of our model. A particular question is 
whether the increase in the virial temperature inside 100 pc 
due to the presence of a SMBH leads to efficient infall to 
smallse radii or not. 

We conclude that a cooling flow can be induced by 
Compton cooling in the UV luminous quasar case, since 
in the clustering scenario of formation of cosmic structure 
quasars are expected to be formed in highly dense massive 
self-gravitating gas clouds (Katz et al. 1994). 
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Figure 1. the Compton cooling time for L\jy = 10 46 erg s 1 and Ltjv = 10 47 erg s 1 and the radiative cooling time for initial gas with 
various temperature; Initial gas distribution is given by equation (12). 



Figure 2. Evolution for model 7.d747. (a), (b), (c) and (d) are the results of distribution of density, velocity, temperature and accretion 
rate, respectively. In (a), the density distribution becomes more steep with time. In the central region, the density increases from 
6 X 10~ 26 gcm~ 3 to 5 X 10 -24 gem" 3 . Since the gas concentrates, the hot gas flows inward from the outer region. As shown in (c), 
in the central region, the temperature decreases from 3 X 10 7 K to 1 X 10 6 K by Compton cooling. Then the flow velocity increases to 
2.4 X 10 7 cms -1 and the accretion rate also increases to 0.4 Mq yr~ 1 in the central region. The gas temperature declines gradually, but 
when the radiative cooling time becomes shorter than the Compton cooling time, the gas temperature declines rapidly. From 10 s yr to 
10 9 yr, the gas flow changes into the quasi- hydrostatic state. 



Figure 3. The same as Fig. 2 , but for model 7.d746 As shown in (a) and (c), at 10 8 yr the density increases and the temperature 
decreases, but the profiles are not more rapid than those in Fig. 2(a) and (c). From 10 8 yr to 10 9 yr, the gas state returns to the initial 
hydrostatic state. For this case, the gas temperature is about 2 X 10 7 K at R c . 



Figure 4. The same as Fig. 2, but for model 5.d747. 



Figure 5. The same as Fig. 2, but for model 5.d7rad. 



Figure 6. The same as Fig. 2, but for model 5.d746. 
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